Topological Hubbard model and its high-temperature quantum Hall effect 
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The quintessential two-dimensional lattice model that describes the competition between the 
kinetic energy of electrons and their short-range repulsive interactions is the repulsive Hubbard 
model. We study a time-reversal symmetric variant of the repulsive Hubbard model defined on a 
planar lattice: Whereas the interaction is unchanged, any fully occupied band supports a quantized 
spin Hall effect. We show that at 1/2 filling of this band, the ground state develops spontaneously 
and simultaneously Ising ferromagnetic long-range order and a quantized charge Hall effect when the 
interaction is sufficiently strong. We ponder on the possible practical applications, beyond metrology, 
that the quantized charge Hall effect might have if it could be realized at high temperatures and 
without external magnetic fields in strongly correlated materials. 



I. INTRODUCTION 

High-temperature superconductivity^ and the quan- 
tum Hall effect (QHEpl have been two of the central 
problems in condensed matter physics of the past three 
decades. The former is related to electrons hopping on a 
two-dimensional (2D) lattice close to (but not at) half 
filling, while the latter focuses on fermions in doped 
semiconductor heterostructures or graphene in a high 
magnetic field. High-temperature superconductors are 
strongly interacting systems, with the potential energy 
about an order of magnitude larger than the kinetic en- 
ergy. In the QHE, the kinetic energy is quenched by the 
external magnetic field. Moreover, interactions are im- 
portant only in understanding the fractional QHE but 
not in understanding the integer QHE (IQHE). 

The possibility that the IQHE could arise in a lattice 
Hamiltonian without the Landau levels induced by a uni- 
form magnetic field was suggested by Haldane in 1988 3 . 
The essence is that, despite the absence of an uniform 
magnetic field, the system still lacks time-reversal sym- 
metry. More recently, it was shown that the fractional 
QHE could also emerge in flat topological bands when 
they are partially fillecP^ (see also ITT)|) . These recent 
developments point to a natural marriage between the 
QHE and strongly correlated lattice systems at a high 
filling fraction. 

In this Letter, we study a quintessential strongly cor- 
related lattice 2D system but with a twist. We consider a 
time-reversal symmetric fermionic Hubbard model in the 
limit of large on site repulsion U compared to the band- 
width W of the hopping dispersion, but with hopping 
terms yielding topologically nontrivial Bloch bands in 
that they each support a quantized spin Hall conductiv- 
ity when fully occupied^. The time-reversal symmetric 
Hubbard model with a single half filled nested Bloch band 
has a charge insulating ground state that supports anti- 
ferromagnetic long-range order^. In contrast, the ground 
state of our time- reversal symmetric Hubbard model with 



topologically nontrivial Bloch bands simultaneously dis- 
plays Ising ferromagnetic long-range order and the IQHE 
at some commensurate filling fraction. The energy scales 
that can be attained in lattice models are typically rather 
high, of the order of atomic magnitudes, i.e., electronvolt. 
If an interacting system with topological bands can be 
found so as to display the IQHE at high temperatures, 
it could be of practical use, as we shall explain after we 
substantiate our claims. 



II. STUDY OF THE TOPOLOGICAL HUBBARD 
MODEL 

We consider spinful electrons hopping on a bipartite 
square lattice A = A\JB with sublattices A and B, where 
each sublattice has N := L, r x L„ sites. The Hubbard 

y 

Hamiltonian with repulsive interactions (U > 0) can be 
written 

H ■= C k n k C k + UY^ n r,t^ n r,i,a- ( 2 - la ) 

fceBZ r a=A,B 

The component cl of the operator- valued spinor ct 
creates an electron with momentum k from the Brillouin 
zone (BZ) of sublattice A and with spin a =f, J,, whose 

Fourier transform c\ a a = N^ 1 / 2 J2keBZ e ~ %k r c k,<r,a is 
exclusively supported on sublattice a = A,B. The 4x4 
Hermitian matrix rl k obeys the time-reversal symmetry 
(TRS) 

U +k = <T 2 H*_ h ° 2 , (2Tb) 

and, owing to a strong intrinsic spin-orbit coupling, the 
residual spin-rotation symmetry (RSRS) 

fh m \ 

H +k = CT 3 ^+fc CT 3 =[ J! ,(4.) > ( 2 - lc ) 
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where the Pauli matrices a ly a 2 , and cr 3 act on the elec- 
tronic spin- 1/2 degrees of freedom. Hence, the two 2x2 
Hermitian matrices h. a ^ with a =^, ^ obey 



AT) 



= h 



-k,Ba> 



Vfe e bz, 



a, /3 = A, B, 

(2.1d) 

because of the condition of TRS (2.1b). Finally, the op- 
erator n r a a = cj, a a c r a a measures the electron density 
on site r in sublattice a and with spin a. 

The Hubbard Hamiltonian defined by Eq. ( |2.1a ) thus 
has a global Z 2 x U(\) symmetry that arises because of 
the TRS ([2Tb} and the RSRS pic]. We are going to 



show that TRS is spontaneously broken while the con- 
tinuous RSRS is shared by the ground state, when this 
Hubbard Hamiltonian acquires suitable topological prop- 
erties. 



It is the choice for the matrix elements h k ^ a/3 



en- 



tering the kinetic energy (2.1c) that endows the Hub- 



bard Hamiltonian (2.1a) with topological attributes. We 
choose 



n k.AB ~ n k.BA - — 



-iw/4 ^ _j_ e +i(k y -k x )\ 
+ e +W4 ( e -ik x +e +ik y j 



At) 

k.AA 



where 



h KBB w k [2i 2 (cos k x - cos k) + 4/i s ] 



:=«£(, + (!-«), kg [0,1], 



(2.2a) 



(2.2b) 



and 



' 1 + cos k x cos k y + [2t 2 (cos k x — cos k x ) + 4/i s ] . 

(2.2c) 

In the noninteracting limit (U = 0), this model fea- 
tures four bands with two distinct two-fold degener- 
ate dispersions ±w k e^. This two-fold degeneracy is a 
consequence of the Kramers degeneracy implied by the 
TRS (2.1b). If we denote the corresponding eigenspinors 
[Xfc a a a) j where A = ±, and choose the normal- 



Xfc.er A 



cr,A' 



= 5 



ization xL ,\X k , 

is diagonalized by using the fcrmionic creation operators 



A, A'' 



Vfc, then the kinetic energy 



fc.(7, A ■ / j *.k,(T,\,a^k,o;a' 
a=A,B 



(2. 2d) 



exactly flat with eigenvalues ±1. The case k = cor- 
responds to a tight-binding model on the square lattice 
that involves only nearest-neighbor (\tA = 1) and next- 
nearest-neighbor hopping (t 2 ) together with a staggered 
chemical potential /i s that breaks the symmetry between 
sublattices A and For k £ (0, 1], longer-range hop- 
ping is introduced. However, we stress that the Hamil- 
tonian remains local for all n e [0, 1] since all correlation 
functions decay exponentially due to the presence of the 
band gajP. 

The topological properties of the lower pair of bands 
are characterized by their spin Chern number 



C s := (C t - CJ /2, 



(2.3a) 



where C a is to be computed from the orbitals of spin-cr 
electrons according to 



C' 



fceBZ 



2tH 



V, A 



(^k,a,— ^/cXfc.cr,— 



(2.3b) 



Time-reversal symmetry implies = — and therefore 
entails a vanishing of the total (charge) Chern number- 
ed := (CL + Ci)/2 of the lower bands. The spin Chern 
number of the lower pair of bands is given by 



a 



i 



AA 



sgnh 



(t) 

(tt,0),AA 



(2.3c) 



Hence, the Bloch bands are topologically trivial when- 
ever | ^2 / A*s I < 1) while the model at half filling exhibits 
the physics of a quantum spin Hall insulator whenever 
1^2/ Ms I 1- I n an °P en geometry, the spin Hall conduc- 
tivity is quantized to the value = eC s / (2ir) where e 
denotes the electric charge of the electron. 

We now consider the system with a repulsive Hubbard 
interaction U > at 1/2 filling of the lower band (1/4 
filling of the lower and upper bands), i.e., with 



xL y = N 



(2.4) 



electrons. In all that follows, we assume that U is much 
smaller than the gap A induced by a strong intrinsic 
spin-orbit coupling between the two pairs of bands. If 
so, we can restrict the iV c -body Hilbert space to the Fock 
space arising from the single-particle Hilbert spaces of the 
lower pair of bands. 

In the limit of flat bands k = 1 and at the commen- 



surate filling fraction (2.4), the kinetic energy (2.2e) at 
fixed spin polarization S := |(cr 3 )| = 0, 2, 
of h/2 has a ground state degeneracy 



./V in units 



as 



Hn 



fcGBZ er= 



t,i A=± 



Ad 



fc.tx.A w k £ k <4v 



(2.2e) 



Hence, the Bloch states created by ct k a x are generically 
spread on both sublattices A and B. We shall consider 
only the case in which these bands are separated by an 
energy gap, i.e., \t 2 \ 7^ |/x s |. The parameter k controls 
the bandwidth of these bands. For k = 1, the bands are 



, N 

■Nes = I JV-|S| 



(2.5) 



The repulsive Hubbard interaction lifts this degeneracy 
whenever any one of these states allows a site of A to be 
doubly occupied with a finite probability. The only two 
states with full spin polarization S = N, 



fceBZ 



.10), 



=t,4, 



(2.6a) 
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are immune to the presence of the Hubbard repulsion. 
More formally, observe that Hamiltonian H — fiN e is a 
positive semidcfinite operator for k = 1, U > 0, and the 
chemical potential /i = — 1. Since 



0, 



=t,4, 



(2.6b) 



the two states (2.6a) belong to the ground state manifold 



of H + N e for any U > 0, t 2 , and /x s . 

We are going to argue that this pair of degenerate Ising 
ferromagnets spans the ground state manifold for any 
U > and \t 2 /n s \ 7^ 1. This is achieved by arguing that 
they are separated from excited states by a many-body 
gap, a departure from the usual ferromagnetism in flat 
bands when full spin-1/2 SU(2) symmetry is not explic- 
itly broken^. First, particle-hole excitations of that 
keep S = N fixed, cost an energy A > and are thus 
gapped. Second, we ask whether excitations of that 
flip one spin (S = N — 2) are gapped as well. Any such 
state can be written as 



V,Q/ = E A k , 4+Q,ff,- < WI*a>. V=tl, 
keBZ 

(2.7a) 

where the center of mass momentum Q is a good quan- 
tum number and thus (<f> a q\§o-,q>) = $q q> if the nor- 
malization X^fcgBZ Afc*A^ = 1 is imposed. One veri- 
fies that (see Appendix) 



(K, Q \(H + N e )\^ <Q 
U 



u- 



N 



E 



V a {q K v 

/ . k A — k—Q,a.~.aAk,cr,— ,a 



fceBZ 



(2.7b) 



where the lowest energy state with one spin flipped is 
characterized by the that minimizes Eq. (2.7b) 



while satisfying the normalization condition. For exam- 
ple, if the single-particle orbitals are fully sublattice po- 
larized, e.g., Xfe CT _ oc (1,0) (topologically trivial), the 

choice = N^ 1 / 2 minimizes Eq. (2.7b) with the 

right-hand side equal to zero. Hence, the fully spin- 
polarized state Y&e) is a gapless ground state in this case. 
On the other hand, let us assume that 

XL,- 9* (M) and <t (0, 1) (2.8) 

holds almost everywhere in the BZ, i.e., up to a set 
of measure zero. In the thermodynamic limit, where 
the sum over k becomes an integral, this delivers from 
Eq. (2.7b) the strict inequality (see Appendix) 



(^ <7>Q \(H + N e )\^ Q )>0. 



(2.9) 



Hence, assumption ( |2.8[ ) is sufficient to show that the 
spin-polarized state is a gapped ground state of the 
Hamiltonian with flat bands in the thermodynamic limit, 
provided that one also assumes that the lowest energy 



states with more than one spin flipped are higher in en- 
ergy than those with one spin flipped. 

Ruling out the possibility that states with many spin 
flips have lower energies than states with few spin flips 
relative to the Ising ferromagnetic ground state is plau- 
sible in the regime when the intrinsic spin-orbit coupling 
generates the largest energy scale (A 3> U). 

Equation (2.8) is a reasonable assumption when the 



Bloch states stem from a band with nonzero (spin) Chern 
number, since the spinor Xk a m &ps out the entire sur- 
face of the unit sphere as k takes values in the BZ. The 
assumption underlying Eq. (2.8) can also be understood 



by constructing the Wannier wavef unctions, centered at 
the lattice point z, of the lowest energy band with spin 
a and Chern number: C„ 



^r, a ,-, a -=i E eih<r ~ Z) X k ,,,-, a - (2.10a) 
feeBZ 

The gauge- invariant part of their spread functional sat- 
isfies 

( to,*,- I r Vo,,,- ) - E K I HVW > I 2 

Z 

> \C a \A c /(2n), 

(2.10b) 

where A c denotes the area of the unit cell. This inequal- 
ity relates the Chern number of the band and the "min- 
imum width" of the Wannier states (see Appendix). In 
particular, in the nontopological phase, one can imagine 
a limit in which the wavcfunction is entirely localized on a 
given sublattice, while the nonzero Chern number in the 
topological phase implies that the Wannier wavefunction 
has amplitudes on both sublattices. 



While Eq. (2.9) is strongly suggestive of the existence 
of a many-body gap A, it does not provide informa- 
tion about its size. To quantify A, we diagonalized the 



model (2.1 ) exactly numerically in the limit of flat bands 
k = 1 at the commensurate filling fraction (2.4). We 
varied the ratio I^/MsI, keeping t| + Ms = 1/2 constant 
to drive the system from the topological to the trivial 
phase. The results are shown in Fig. [T] First, they sup- 
port the assumption that all states with more than one 
spin flipped are higher in energy than the many-body 
one-spin-flipped gap provided \t 2 /n s \ > 1. Second, we 
find A s» 0.3U as an extrapolation to the thermody- 
namic limit for /i s = deep in the topological phase, 
while A monotonically decreases toward a much smaller 
nonvanishing value for t 2 = in the topologically triv- 
ial phase set by the unit of energy = 1. Finally, 
it should be noted that neglecting the states from the 
upper band of the noninteracting Hamiltonian delivers 
the correct excitation many-body gap A not only in the 
aforementioned limit U <S A , but also under the weaker 
condition A < A , if the limit of flat bands is taken. 

Deep in the topologically nontrivial regime l^/Msl 
1, the states and are degenerate ground states 
related by TRS for any finite N. They are separated 
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FIG. 1. Numerical exact diagonalization results for flat 
bands k = 1 at the commensurate filling fraction (2.4 1. 



Markers show the energy of the lowest state in different sec- 
tors of total spin S (in units of h/2) measured with respect 
to the ground state energy for L x = 3, L y = 4. Here, 
g :— (2/7r)arctan|/i s /t 2 | so that g > 0.5 and g < 0.5 cor- 
respond to the trivial and topological single-particle bands, 
respectively. Since there is only one state in the fully po- 
larized sector | ,S| = 12, the difference between the asterisks 
and the squares is the many-body excitation gap A(g). The 
thick blue line shows the extrapolation of A(g) to the ther- 
modynamic limit. In the inset, exact diagonalization in the 
sector with one spin flipped away from the fully polarized sec- 
tor is presented for p s = 0, t 2 = l/\/2 and L x = L y ranging 
from 6 to 30. The straight lines are a guide to the eye and 
make evident an even-odd effect in L x = L y . Deep in the 
topologically nontrivial regime g <C 0.5, we observe a sizable 
A(g <C 0.5). The topologically trivial regime g > 0.5 is also 
characterized by a gap A(g > 0.5) in the sector with one spin 
flipped away from the fully polarized sector, however this gap 
is much smaller than A(g <C 0.5). We refer the reader to the 
Appendix for a discussion of the regime A(g > 0.5). 



from their excitations by a gap that survives the thermo- 
dynamic limit N — >■ co. Spontaneous breaking of TRS 
takes place in the thermodynamic limit N — ¥ oo by se- 
lecting the ground state to be |^), say. It is then mean- 
ingful to discuss the quantized electromagnetic response 
of |^*), since TRS is spontaneously broken. The trans- 
verse charge response a^ y of |^) is proportional to the 
many-body Chern number C,^ >. The latter takes into 

account the occupation of the Bloch states 4 . Since all 
Bloch states of the lower band with spin a are occupied 
in 1^^.), while all Bloch states with spin J, are empty, 
C^j, \ = Cj.- Hence, the ground state has the quantized 
Hall response 



|(7» | = |C f | x e 2 /h = e 2 /h. 



(2.11) 



Remarkably, the selection by the repulsive Hubbard 
interaction of a ground state supporting simultaneously 
Ising ferromagnetism and the IQHE is robust to a sizable 
bandwidth as is suggested by numerical exact diagonal- 
ization. As shown in Fig. ([2]), the fully spin-polarized 
state \^ a ) remains the gapped ground state of the sys- 
tem up to a bandwidth W/U « 0.7. 
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FIG. 2. (Color online) Numerical exact dia gonalization re- 
sults at the commensurate filling fraction \2A\ as a function of 
the bandwidth W for L x = 3, L y — 4. Plotted is the energy 
of the lowest state in different sectors of total spin S (in units 
of h/2) measured with respect to the ground state energy in 
the topological phase with p s — 0, t 2 = \/\f2. The ground 
state is gapped and fully spin-polarized for W/U < 0.7, while 
it is unpolarized for W/U > 0.7 (see Appendix). 



III. PRACTICAL APPLICATIONS 

So what is it good for, a material with a QHE at room 
temperature wi thou t applied external magnetic fields be- 
sides metrolog y* 14 * 15 *? First, we recall that the quantiza- 
tion of the Hall resistance and the accompanying van- 
ishing of the longitudinal resistance is exact only at zero 
temperature. The longitudinal resistance increases expo- 
nentially fast with increasing temperature^. However, if 
a QHE with gaps of the order of hundreds of meV or even 
eV scales could arise in a strongly correlated lattice ma- 
terial, exceptionally low resistivities could be attained. 
The resistance of a Hall bar depends on its aspect ratio 
and the Hall angle 5 = &ictan{p xy / p xx ^^, but for long 
systems ("wires") near the quantized regime, the lon- 
gitudinal resistance scales as R xx = L/W p x , r , and the 



2D resistivity p x 



K 



A / T where A is the exci- 



tation gap. For gaps of the order of 100 meV to 1 eV, 
one would obtain room temperature 2D resistivities from 
Pxx ~ 10 3 to p xx ~ 10~ 13 f2, respectively. Obviously 
the exponential behavior is responsible for this gigan- 
tic range. Small as they are, these are not perfect con- 
ductors. For a benchmark, we consider the conductivity 
of copper at room temperature per atomic layer. Us- 
ing the value for the 3D resistivity of copper at 20° C of 
Pol = 1.68 x 10~ 8 fim^and that the lattice parameter 
for FCC lattice is 3.6lA, we obtain = 93.3 fl. There- 
fore, for gaps above A w .2eV, the Hall system starts 
to be better conducting than copper at room tempera- 
ture, and for A ss .3eV it is already almost 3 orders of 
magnitude better conducting than copper. 

The RSRS (2.1c I is not exact in practice. For example, 
a Rashba spin-orbit coupling violates this RSRS. How- 
ever, our analysis of transport at room temperature still 
applies provided the characteristic energy scale associ- 
ated to the breaking of the RSRS (2.1c) is much smaller 



than the largest energy scale A induced by the intrinsic 
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spin-orbit coupling. Mater ials that realize a 2D Z 2 topo- 
logical band insulato r 18 * 19 ^ with a band gap A are thus 
candidates to realize a QHE at room temperature if (i) 
the band gap is larger than the correlation energy and 

(ii) the chemical potential can be tuned to half-filling of 

(iii) a reasonably flat valence band. HgTe quantum wells 
with an inverted band structure realize a 2D Z 2 topo- 
logical band insulator with a small Rashba coupling 20 21 . 
The design of a material with the functionalities (i)-(iii) 
has been proposed in Re f. [H . Cold atoms trapped in an 
optical honeycomb latticd^lMI might offer an alternative 
to realizing the topological Hubbard model discussed in 
this Letter. 

We would like to close by mentioning that examples 
such as the topological Hubbard model discussed in this 
Letter, as well as lattice models displaying the fractional 
quantum Hall effect studied in Refs. EHI1 could serve as 
benchmarks for numerical methods of fermionic models 
in 2D such as dynamical mean-field theory and meth- 
ods based on tensor product states^. In contrast to 
the single-band repulsive Hubbard model, for which lit- 
tle is known exactly at fractional filling, the topological 
Hubbard model (2.1 ), because of the nonvanishing Chern 



numbers of its bands, leads to much better understood 
(topological) ground states. It can thus serve as a yard- 
stick for the performance of these methods. 

This work was supported in part by DOE Grant No. 
DEFG02-06ER46316 and by the Swiss National Science 
Foundation. 



in the N e -many-body state 



\**,q) = E 4 Q) 4 +Q , ff ,-4, CT ,-I* CT >, o =t4, 

feeBZ 

(Alb) 

which has one spin flipped as compared to the - up to 
time-reversal symmetry and the global gauge phase fac- 
tor - unique normalized N e -many-body state with full 
spin polarization 



(Ale) 



fceBZ 



Here, Q denotes the center of mass momentum and the 
orthonormalization condition 



(*<t,qI*<t,Q'> =<W 



(Aid) 



Appendix A: Intermediary steps 
1. Derivation of Eq. fiZTFbty 



enforces the normalization 



For the limit of flat bands k = 1, we are going to 
compute the expectation value 



(^ Q \(H + N e )\$„ tQ ) 



(Ala) 



E 4 Q) *4 Q) = i- 



fceBZ 



(Ale) 



We rewrite |3> ct ,q) as 



1$ 



<T,Q> 



E E X k+Q,S -,a C k- 



fceBZ 



--A,B 



^ A iQ) 

E k 



E 



0=A,B 



E E 

a,/3=A,B r,r'eA LfcGBZ 

E E <2„< 

a,P=A,B r,r'eA 



N 



* _ -i(k+Q)-r +ik-r' 

Afc+Q,er, — ,aAfc,£r,— ,a° ' 



(A2a) 



/3 ^r,a,a^r' ,a,f3 I *ct/j 



where we have introduced the short-hand notation 



M, 



A 



(Q) 



r,a,r',/3 " ^ j\j 

fceBZ 



E^k_ * -i{k+Q)-r +ik-r' 



(A2b) 



G 



in terms of which 



feA S=A,B 



= u E E \KU 

r,r,r'eA a,f3,S=A,B 



r,r,r'eA a,/3,S=A,B 

u E E HU,\ 2 -uY E \ M $lr, a 

r,r' GA a,ft=A,B rGAa=A,B 



(A3) 



Using the definition ( A2b I , the first term becomes 



^(Q)^(Q)* 



V \ " kf(t) V V \ " ^k * * i(fc _ fc ').(y_ r ') 

/ , / , \ 1U r,a,r',l3 ~ Z-^ 2^1 JJ2 \ A: ■() .. . ■ U.\ ... U • Q i\k\\ 



r,r'eA a,f3=A,B 



r,r'GA a,p=A,B fc.fc'GBZ 



^(Q)^(Q)* 



EV ^ "^fe ^fe' * * 7V7"2 C 

/ . Xfc+Q,4.,-,/3Xfe,t,-,aXfc'+QJ.,-,^Xfc',t,-ia-^' 
a,/3=A,B fc,fc'eBZ 



V A {Q 'A iQ] * V v* v 

A k A k Xfc+Q,4.,-,/3X. 



fc,GBZ 

^ E 4 Q) 4 Q> 

fc,GBZ 
1, 



k Q \ J j E Xfe,t,-,o;Xfe,t,-,C« 

i a—A,B 



(A4) 



where both the normalization conditions on -A^ i n Eq. (Ale) and on x& o- _ a have been used in the last two lines. 
The second term in Eq. ( A3 ) becomes 



E E K 

reA a=A,B 



fit) 

a.r.a 



A (Q) A {Q)* 



E E E - TV 2 Xk+Q,i,- 



r£A a=A,B k,k'£BZ 



tXk,-l,-,aXk'+Q,l,-,aXk',t,- 



E E 



A k A k' „* 



a=A,B k,k'£BZ 



N 



<}... .„\k \ .,,\k -(j... ,,\k\\ .,, 



- E 

a=A, 

- y 



EAQ) * 
A k Xk+Q,i,-,aXk 



t,-,a 



feeBZ 



(A5) 



=.4.B 



/ / A k X-h-Q,t,-,aXk,t,-,a 



keBZ 



where we used the identity Xk a \ a ~ X-k s x a which follows from the time-reversal symmetry of the Hamiltonian. 
Putting everything together, we obtain Eq. (2.7b I of the Letter, 



($ a , Q \(H + N e )\$ atQ ) = U-^ ]T 



N 



-A.B 



V a {q K v 

/ j h A -fc-Q,(7,-.a A fc , cr, — , a 



(A6) 



2. Derivation of Eq. (2.9) 



Before starting with the derivation of the inequality 
Eq. (2.9), let us establish the following inequality 



for a, b, c, d £ [0, 1]. In particular, the equality only holds 
when either a = b = c = d = or a = b = c = d= I. 
To show this, one can rewrite the inequality in terms of 



1 > abed + x/l-a 2 \/l-b 2 (A7) 
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trigonometric functions with angles a, (0,7, 8 € [0, 7r/2] 



strict inequality 



1 > sin a sin /3 sin 7 sin S + cos a cos f3 cos 7 cos 5 

— i [cos(a + /3) cos(7 + (5) + cos(a — /?) cos(7 — 5)] . 
1 (A8) 

Hence the inequality holds and the equality is true only 
if either a = /3 = 7 = (5 = 0ora = /3 = 7 = (5 = 7r/2as 
announced. 

We will now show that in the thermodynamic limit the 



iV 



V a {q K v 

/ ^ k A-fc-Q,(T,-,QAfc,ff,-,Q 



fcSBZ 



holds under the assumption that Eq. (2.8) is satisfied 



> 



(A9) 



almost everywhere in the BZ, i.e., up to a set of measure 
zero. 



We rewrite, for very large N, 
U 



N 



E 

a=A,B 



V a {q K v 

/ j k A- — k — Q,cr,— ,aA-k,(T,— , 



feGBZ 



NU 

< NU 

< NU 

< NU 
= U. 



d 2 k 



d 2 k' 



(Q) a(Q)* 



Ar a 



(2tt) 2 J (2tt) 2 fc k 



^ ^ X— fc— Q,a, — ,a 



Xfc.cr, — ,aX— fc' — Q,tT, — ,aXfc' ; CJ", — ,Q 



d 2 fc 



d 2 fc' 



(2tt) 2 y (2tt) 



.4 



(Q) 



.4 



a=A,B 

(Q)* 
fe' 



E 



I X— fc— Q,(7,— ,q I |Xfe,o\ — ,ct I I X— fc' — Q,tx,— ,a I IXfc'.tr,— ,a I 



(A10) 



d 2 fc 



d 2 k' 



(2tt) 2 y (2vr) 
cPfc (Q) 



.4 



(Q) 



(Q)* 



4? 



(2^) 2 



Due to the normalization of the eigenspinors, the assumption (2.8), and the inequality (A7), 

| X— fc— Q.cj, — ,q I |Xfe,iT,— 



E 



|X*fe'-Q,cr-,a| |XfcV,-,a| < 1 ( A11 ) 

holds for almost all k, k' . This allows the use of the strict inequality in the third-last inequality from the chain of 
inequalities (A10). 



r 



The penultimate inequality from the chain of inequal- is non- vanishing under the assumption (2.8). 
ities (A10) is a consequence of Holder's inequality 



d 2 k 



A 



(Q) 



< 



d 2 k 4Q) 2Xl2 



(2tt)' 



(A12) 



The last inequality in the chain of inequalities ( A10) fol- 
lows from the representation 



d 2 k 



A 



(Q) 



1 

N 



3. Proof of Eq. ( |2.10b[ ) 

We are now going to work exclusively with the single- 
particle eigenstates of the band Hamiltonian (2.2c). 
We start with the completeness relation 



(A13) 



of the normalization ( |Ale ) in the thermodynamic limit 
In summary, we have shown that 



E E \r,a)(a,r\ = l 



(Al6a) 



r£A a=A,B 



(K.q\(H + N c )\^^ q ) >0 

holds in the thermodynamic limit under the assump- 
tion (2.8). In other words, the many-body gap 

A a Q := ($ CT|Q \H\ $ CT , Q ) - \H\ * a ) > (A15) 



on the bipartite lattice A = A\JB. The normalized Bloch 
(A14) states \tp k are defined in terms of the single-particle 
eigenstates Xk a \ a by their overlaps 



-i k r 



Xfc,er,A,. 



(Al6b) 



for any given (Ising) spin a =t,i and band A = ±. The 
normalization and phase factors are chosen so that the 
orthonormality condition 



(<Ph, t r,x\'Ph,a,X') = S X,X' 



(Al6c) 



with A, A' = ± holds for any given wavenumber k G BZ 
and any given spin a =t,|- The overlap (Al6b) is in- 



variant under the translation k > k + Q for any Q that 
belongs to the reciprocal lattice of sublattice A owing to 
the periodicity 



Xfc,cr,A,a Xfc+Q,o\A,. 



(Al6d) 



Wannier states are defined in terms of the Bloch states 
by the unitary transformation 



1 



E 

fceBZ 



e _ifc -*l¥W) 



(A17a) 



for any given unit cell of sublattice A labeled by z and 
any given spin a =T,4- The orthonormality ( Al6c[ ) of 
the Bloch states thus carries over to the orthonormality 



l1>z,<r,\\i>z,<r,X') = S X,X> 



(Al7b) 



of the Wannier states for any given unit cell labeled by 
z of sublattice A and any given spin <r =t,|- For any 
A and a, the representation of the Wannier state on the 
bipartite lattice A = A U B is the overlap 



1p r 



,x,a (r,a\ip Xt<TiX 
1 

~~ y/N 



e-* kz (r,a\<p kt<T>x ) 



keBZ 



- E 



+ih-(r-z) 
c Afc.er.A,. 



(A17c) 



fceBZ 



for any given unit cell of sublattice A labeled by z and 
any given spin a =f, ^. The set of Wannier spinors 



Vv, 



z,a. A 



r,z,a,\,a ) 



(Al7d) 



resolves the identity since 



E< 

reA 



A^V,*',ct.A 



A, A' 



(A17e) 



for any given spin a =f, \.. We want to estimate the 



profile in space of the Wannier states (Al7a 



For that purpose, we consider the spread functional 



R*l :=(^=o )CT ,Ak 2 |^=o, CT ,A> 



(A18) 



Here, we are assuming, for simplicity, that there is only 
one band A that is occupied. If more bands are occupied, 
we have to carry out a summation over all the occupied 

(21 

bands. Observe that by translational invariance R} { is 
left unchanged under the global translation r, z —> r + 
R, z+ R for any lattice vector R. Hence, the choice z = 
in Eq. (A18) can be done without loss of generality. 



We rewrite Eq. (A18), following Ref. [13], as 



R 



(2) 
cr,X 



p(2) 

K l\<T,X 



m -.= 



I|<T,A (7, A' 

(i } Z=0,a,x\ r2 \'>Pz=0,a,x) 

E i (^=0,^1 

z' 

E \ ^z=0, a ,x\ r \^z-,a,X 



(Al9a) 

| 2 , (A19b) 
; . (Al9c) 



(2) ~ (2) 

Both i?V A and R a x are non-negative quantities that can 
be expressed in terms of fc-space summations as follows. 
First, 



^I|2,A - E AT E [(^j^Lo-.a) {pk } Xk,a,xj {pk^k,a,\ Xk,a,x) (xk,a,\ ^Xfe.ov 
j=x,y keBZ 

~ E AT E [( ®kXk,cr.\ I ^ Xfe.tr, a) — ( dkXk,o,\ I Xfc.cr, x) ( Xh,a,X I ^Xfe.er.A ) 
j=x,y keBZ 

= E ^7 E^^-^aI ( 1_ I Xfe.^A M Xfe.o-.Al) l 5 fe.Afe, CT ,A> 
j=x,y keBZ 



E E ^fe^fe.o-.AlQfe.a.Al^fe.Xfe.a.A) 

3=a,j/ fceBZ 

4 E tr [Sfc,«r,A] ' 



(A20a) 



fceBZ 
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where 



and 



Qk : a,X '■— 1 I Xk,a,\ ) (Xfe,<7,Al 



(A20b) 



is the single-particle projector operator on all the Bloch 
states orthogonal to | \k a x ) > tne 2x2 matrix [g k a x ] 
has the components 



9k,a,X\i,j : ~ R- e (^fe.Xfe,<7,A I Qk,cr,x\ ^fc.Xfe,er,A/J ' 

(A20c) 

labeled by the Euclidean indices i,j € {x,y} of two- 
dimensional space, and tr denotes the trace over i,j £ 
{x,y}. Second, 



R 



(2) = \: 
<r,A L^i ]\r 

j=x,y k 



V ( A k,<T. 



A|J A cr,X\j 



(A21a) 



where 



.4 



k,a,X\j : ~ * \Xh,a,\\ ®kjXk,a,x)> 



A 



k 



(A21c) 



denote the Berry connection and the average of the Berry 
connection, respectively. Under the gauge transforma- 
tion 



IXfc,<T,A / 



IXfc,cr,A /) 



(A22) 



(2) ... . ~ (2) 

it can be shown that R), , remains invariant while R„ { 



does not. 

We will now establish a lower bound on the gauge in- 

(2) 

variant quantity R\^ x - To this end, we notice that, since 
Q k a x is a projection operator with eigenvalues 0, 1, then, 
for any single-particle state a A ), it follows that 



< "I'a-.^.a Qk,?.x 'I'a- , A 
In particular, if we choose 



I* 



±, fc.tr, A , 



|oLXfc,a,A) ±i \ d k„Xk,a. 



A/> 



(A21b) the inequality (A23) delivers 



(A23) 
(A24) 



tr [9k,<r,\] > Ti (( d K Xk,a,x I \Xk, a ,X ) - ( d k y Xk,*,X I d k m Xk 



a.X 



Summing (A25) over k and taking the thermodynamic limit N — > oc, delivers 

A,. 



R (2) 



k,A ( 2?r )2 



BZ 



> ± 
= ± 



A„ 



d 2 k 



2n 7bz 27ri 



2tt 



A ' 



, dk x Xk,a,X I ®k Xk,a,x) ~ ( ^fe Xk,a.X I Xfc : cr,> 



(A25) 



(A26) 



where A c is the area of the unit cell of sublattice A and 
C a x are the spin- and band-resolved Chern numbers. 

(2) 

bmcc R\V , > 0, it then follows 

1|£T,A — 1 



R 



(2) 



A 

> ^ IC 



o-.Ah 



(A27) 



as a lower bound, proportional to the band Chern num- 
ber, on the gauge invariant part of spread of the Wannier 
states ( A17a). It was shown in Ref. [fT5] that, in the topo- 



logical phase of the Haldane modeP, R^ a x is finite while 

R x displays a logarithmic divergence which is related 
to the non-zero Chern number of the bands. 



Appendix B: Competing Phases in the Hubbard 
model 



Any translationally invariant electronic single-particle 
Hamiltonian is described in momentum space by the dis- 
persions and Bloch eigenfunctions of its bands. Given 
the hopping amplitudes of a tight-binding model in posi- 
tion space, the dispersions and Bloch eigenfunctions are 
uniquely determined. Once interactions are included, the 
many-body ground state of the model depends crucially 
on the properties of the dispersions (bandwidths; shapes 
of Fermi surface, e.g., nesting; ...) and the properties 
of the Bloch states or Wannier states (Chern numbers, 
localization in position space, ...). The power of our ap- 
proach, the flattening of the band with the factor w k 
from Eq. (2.2), is to disentangle the properties of the dis- 
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FIG. 3. (Color online) Comparison of the topological (top; 
t 2 — l/v2> /J. s = 0) and non-topological (bottom; t 2 = 0, jj, s = 
l/v2) single-particle model, (a) The shaded area represents 
the Fermi see of the lower band at the commensurate filling 
fraction N e — N when k < 1. (b) The eigenspinor Xfc,o-> when 
interpreted as a point on the surface of the unit sphere, swipes 
out the full surface of this sphere (a small portion of this 
sphere near one pole) as fc takes values everywhere in the BZ 
for the topological (non-topological) band structure, (c) The 

in real space indicates 



spread of the Wannier states (A17c 



their delocalized (localized) character for the topological (non- 
topological) band structure. 



persions from those of the spinors to some extend and to 
study their effects individually. 



Possible phases 
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FIG. 4. (Color online) Numerical exact diagonalization re- 
sults for flat bands k = 1 for the system size L x — L y = 25 
at the commensurate filling fraction N e = TV. Plotted is the 
value of the many-body excitation gap A/U between the en- 
ergy of the fully spin polarized state S = N and the lowest 
energy state in the sector with one spin flipped S — N — 2 as a 
function of g := (2/7r)arctan|pi s /f 2 | and the nearest-neighbor 
hopping parameter t 1 . The thick red line corresponds to the 
parameter choice made for the plot in Fig. 1 of the Letter. As 
discussed in the text, the spinors of the Bloch band become 
fully sublattice polarized in the limit t\ — > and the fer- 
romagnetic ground state becomes degenerate in energy with 
states from other spin sectors (A — >■ 0). 



Let us list the anticipated ground states for some lim- 
iting properties of the spinors and the dispersions for 
N e = N electrons when the two pairs of bands are sepa- 
rated by a gap in the absence of interactions. 

1. (Almost) flat band and sublattice-polarized spinors. 
If the band is completely flat (as = 1), the ground 
state is macroscopically degenerate in the absence 
of interactions. Moreover, when the Bloch states 
are fully sublattice polarized, i.e., when x\ — (1) 0) 
say, then the kinetic energy becomes diagonal in 
position space [Vvzo-- = ^rz (M)]- Hence, any 
state with no less and no more than one electron 
per site on sublattice A will be a ground state of the 
Hamiltonian, regardless of the spin-configuration 
of the electrons. Therefore, the ground state re- 
mains macroscopically degenerate, in spite of the 
presence of the Hubbard interaction. This degen- 
eracy can be lifted by introducing a small but finite 
bandwidth, while keeping the full sublattice polar- 
ization of the spinors intact. Then, the low-energy 
degrees of freedom of the model map on the con- 
ventional one-band Hubbard model at half filling 



with t -C U, i.e., on the t-J model. An antiferro- 
magnetic ground state will be selected. An alter- 
native way to lift the macroscopic degeneracy is to 
add longer range interactions, in which case other 
ground states may be stabilized^. Note that the 
fully sublattice polarized spinors discussed here im- 
ply that the non-interacting model is topologically 
trivial. 

Flat band and not sublattice polarized spinors. As 
was shown in the Letter, if the spinors are suffi- 
ciently spread between the two sublattices, i.e., if 
condition (2.8) is met, the ground state in the limit 



of flat bands is an Ising ferromagnet with full spin 
polarization. 

3. Strongly dispersing unnested band. For a strongly 
dispersing band without a nested Fermi surface, the 
Fermi liquid ground state is expected to be stable 
against sufficiently small repulsive Hubbard inter- 
actions. One generic instability at finite interac- 
tion strength is Stoner ferromagnetism which, in 
contrast to the flat band ferromagnetism, does not 
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FIG. 5. Numerical exact diagonalization results at the com- 
mensurate filling fraction N e = N as a function of the band- 
width W for = 3, L y — 4. Plotted is the energy of the 
lowest state in different sectors of total spin S (in units of h/2) 
measured with respect to the ground state energy, (a) Topo- 
logical phase with fj, a = 0, t 2 = l/\/2. The ground state is 
gaped and fully spin-polarized for W/U < 0.7, while it is un- 
polarized for W/U > 0.7. (b) Topological trivial phase with 
H s = l/v2, t 2 = 0. The unpolarized ground state appears 
already for very small values of W/U. 



automatically yield full spin polarization. 

Strongly dispersing nested band. If the Fermi sur- 
face is nested, arbitrary small Hubbard interactions 
destabilize the Fermi liquid ground state toward an 
antiferromagnet . 



2. Realization in our model 

The shape of the Wannier states is controlled by the 
parameters t 2 and /i s in our model. If they are topolog- 
ically non-trivial (\t 2 //J, s \ > 1), the Wannier states are 
relatively wide spread over several lattice sites as illus- 
trated in Fig. [3] (c) and by Eq. (A27|. In contrast, deep 
in the topological trivial regime \t 2 /n s \ <C 1, the Wannier 
states become more and more localized. Intuitively, their 
localization is more pronounced the smaller the area cov- 
ered by the spinor Xk a x on the un ^ sphere as k takes 
values everywhere in the BZ [Fig. [3] (b)]. 

In our model, the full sublattice polarization is 
achieved by taking the limit t x — > while \t 2 /fj, s \ < 1 
(this limit is not shown in the Figures of the Letter). In 
the case of flat bands, this limit is the transition from 
phase [2] to phase [T] from the list above. As shown in 
Fig. gj the many-body gap A of the flatband ferromag- 
netism collapses regardless of the value taken by the ratio 
\t 2 / n s \ as t x — > 0. The lowest energy states of each spin 
sector become degenerate in this limit as the full SU(2) 
spin- 1/2 symmetry of the Hubbard interaction is recov- 
ered in the limit t x — > 0. 

We studied the effect of finite band width both in the 
topological trivial (|£ 2 /m s I < 1) an d non- trivial sector 
(\t 2 / n s \ > 1) sector (Fig. [5j . In both cases, an interpo- 
lation between the phases |2] and [4] follows from changing 
the ratio W/U [Fig. |3ja) illustrates the nesting of the 
Fermi surfaces for the parameters chosen.]. The differ- 
ence between topological non-trivial and trivial case is 
quantitative: While the many-body gap A protects the 
Ising ferromagnetic state against the effect of a sizable 
band width in the former case, the transition to the an- 
tiferromagnetic state occurs already for very small W/U 
in the latter. 
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